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CONCEPT  EVALUATION  FOR  OCEAN  SHEAR  PROFILER 


Introduction 

There  presently  exist  several  R&D  devices  that  provide  ocean  shear  profile  information, 
either  in  a  moored  or  a  free-fall  mode  of  operation.  It  would  be  exceedingly  useful  to  have  a 
device  that  could  provide  shear  profile  information  in  the  over-the-side  mode  commonly  used 
for  acquiring  ocean  station  hydrographic  data.  Present  instruments  that  are  operated  in  this 
mode  suffer  severe  degradation  of  performance  due  to  ship  motions.  The  purpose  of  this  note 
is  to  describe  a  new  concept  that  should  be  useful  for  this  purpose  but  which  still  requires  a 
thorough  engineering  design  study  and  prototype  construction,  test,  and  evaluation  to  be  certain 
that  instrument  noise  criteria  can  be  met. 

Present  Technology 

The  spatial  and  temporal  variability  in  upper  ocean  physical,  chemical,  and  biological 
parameters  is  governed  by  dynamical  processes  that  occur  on  many  space  and  time  scales.  The 
unraveling  of  these  processes  requires  detailed  information  on  the  fluxes  of  kinetic  and  poten¬ 
tial  energy  in  the  upper  ocean  layers.  To  meet  these  information  requirements,  new  technology 
has  provided  the  oceanographer  with  sophisticated  instruments  that  measure  water  temperature, 
conductivity  and  velocity,  as  well  as  other  less  dynamically  important  parameters.  In  the  last 
decade,  improvements  in  instrumentation  have  provided  the  CTD  profilers  that  presently  yield 
accurate  conductivity  and  temperature  profiles  on  a  survey  basis.  However,  the  accurate  meas¬ 
urement  of  water  velocity  profiles  is  only  presently  being  attempted  with  specialized  research 
Manuscript  submitted  December  10,  1979. 
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instruments.  There  are  several  types  of  these  profilers  that  have  shown  promise,  but  they  are 
not  in  routine  use  on  survey  vessels.  They  can  roughly  be  grouped  into  mooring,  free-fall  and 
underway  types. 

The  tethered  cyclesonde  of  VanLeer  et  al  (1974)  is  a  current  meter  that  has  a  controllable 
buoyancy  element  and  is  attached  by  a  sliding  connector  to  a  mooring  cable.  At  a  controllable 
rate,  it  slides  up  and  down  the  cable,  yielding  repeated  profiles  at  the  mooring  location.  How¬ 
ever,  the  mooring  float  is  constrained  to  be  well  below  surface  effects  (cf.  Gould,  Schmitz  and 
Wunsch,  1974)  to  prevent  contamination  of  the  velocity  records  by  mooring  motions,  so 
uncontaminated  near-surface  profiles  are  not  attainable. 

The  free-fall  types  are  launched  from  a  research  vessel,  they  descend  to  a  prescribed 
depth  or  the  bottom,  ascend  again,  and  are  recovered.  The  purpose  of  having  them  untethered 
is  to  prevent  contamination  by  ship  or  mooring  motions,  but  of  course  station  time  and  ease  of 
handling  suffer  appreciably.  There  have  been  several  types  of  profilers  in  this  class.  Sanford 
(1975)  uses  an  EM  device  that  measures  the  voltages  induced  by  the  movement  of  the  sea  and 
the  instrument  through  the  geomagnetic  field.  Rossby  (1969)  has  acoustically  tracked  a  freely 
falling  dropsonde,  and,  as  with  Sanford's  instrument,  the  vertical  resolution  is  the  order  of  10 
meters.  Simpson  (1975)  has  used  a  tiltmeter  to  deduce  smaller  scales  down  to  order  30  cm, 
and  Osborne  (1974)  has  obtained  vertical  resolutions  to  1  cm  by  measuring  the  force  on  a  spin¬ 
ning  cantilevered  vane.  These  increased  resolutions  come  only  at  the  expense  of  information 
on  the  longer  scales.  For  ease  in  handling,  Sanford  has  been  developing  an  expendable  profiler 
which  is  based  upon  a  combination  of  the  technology  in  his  EM  profiler  and  in  the  XBT  drop 
body  and  wire  link.  This  of  course  has  the  advantage  of  being  deployable  while  underway. 
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Finally,  underway  profilers  using  acoustic  Doppler  backscatterometry  are  being  developed 
which  are  modifications  of  commercially  available  bottom  speed  indicators  for  supertanker  navi¬ 
gation.  These  profilers  may  be  useful,  but  they  potentially  are  very  sensitive  to  surface  wave 
induced  ship  motions  and  to  the  density  of  volume  scatterers  in  the  water  column.  They  also 
practically  are  limited  to  the  top  few  tens  of  meters  of  water,  and  quite  possibly  will  not  be  use¬ 
ful  in  the  present  context. 

There  is  a  need  for  a  profiler  that  can  be  deployed  in  the  usual  manner  in  which  the  CTD 
is  deployed  and  which  has  comparable  resolution.  In  fact,  it  should  use  the  same  hydrographic 
cable  so  that  both  the  density  and  velocity  profiles  can  be  acquired  at  the  same  time.  It  is  the 
purpose  of  this  note  to  identify  what  appears  to  be  a  workable  system  that  is  relatively  immune 
to  ship  motions. 

Ship  Motions  and  Current  Meters 

The  problem  with  deploying  traditional  current  meters  on  hydrographic  wire  to  obtain 
velocity  profiles  is  sensor  motion  due  to  coupling  with  the  ship.  There  are  two  types  of  motion 
that  cause  problems.  First,  any  vertical  motion  such  as  raising,  lowering,  or  ship's  roll  causes  a 
vertical  component  of  flow  past  the  meter.  In  mechanical  meters,  this  "pumps  up"  the  energy 
level  because  of  nonlinear  coupling  of  the  longitudinal,  up/down  motions  to  the  near-horizontal 
motions  in  the  ocean.  This  problem,  in  principle,  could  be  eliminated  by  nonmechanical  meas¬ 
urement  systems  such  as  the  Neil  Brown  acoustic  current  meter  (N.B.I.S.,  1978).  Second,  drift 
of  the  ship  drags  the  instrument  horizontally  through  the  water  at  a  rate  that  can  be  much 
greater  than  the  signal  being  measured.  This  causes  an  indicated  velocity  that  may  bear  l.ule 
resemblance  to  the  real  one,  and  it  also  tilts  the  instrument  so  that  only  a  component  of  the 
horizontal  velocity  is  registered. 
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The  effect  of  both  of  these  motions  can  be  eliminated  by  a  judicious  choice  of  the  param¬ 
eter  being  measured.  In  this  instance,  we  will  argue  that  the  direct  measurement  of  fluid  vorti- 
city  yields  a  useful  technique  for  measuring  the  required  water  velocity  information  while  at  the 
same  time  suppressing  the  noise  caused  by  instrument  motion. 


Theory 


The  instrument  should  be  designed  to  measure  the  horizontal  components  of  vorticity 

.  Bu_  _  Biv 
tjr“  Bz  Bx 


r  _  Bv  Bw 
y  Bz  6y 


(1) 


where  (x,y)  are  horizontal  axes,  z  is  the  vertical  axis,  and  (u,v,w)  are  the  two  horizontal  and 
the  vertical  component  of  velocity.  The  meter  will  be  concerned  primarily  with  horizontal  velo¬ 
cities,  since  the  vertical  velocity  is  appreciable  only  in  small  scale  motions.  In  fact,  an  approxi¬ 
mation  in  which  the  vertical  motions  are  entirely  neglected  is  well  grounded  in  the  use  of  all 
previous  current  meters.  However,  since  vorticity  is  to  be  measured,  the  relation  between  the 
vorticity  and  the  vertical  profile  of  horizontal  current  is  worth  reconsidering. 


Low  frequency  (mesoscale  and  below)  motions  do  not  have  appreciable  vertical  com¬ 
ponents,  but  the  higher  frequency  motions  may.  In  the  important  internal  wave  band,  the  vert¬ 
ical  velocities  can  be  appreciable.  By  continuity,  the  vertical  velocity  is  related  to  the  horizontal 
one  by 


w—kJkjU,  (2) 

where  k„  and  k:  are  the  horizontal  and  vertical  wave  numbers  of  a  plane  wave.  By  the  disper¬ 
sion  relation, 


(3) 
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where  N  is  the  local  buoyancy  or  Vaisala  frequency  and  /is  the  inertial  frequency,  so  that 


w/u 


G>2-/2 


1/2 


(4) 


Thus,  we  conclude  that  waves  with  frequencies  near  th  e  local  Vaisala  frequency  have  vertical 
velocities  as  large  as  or  larger  than  the  horizontal  velocities.  However,  there  is  not  much 
energy  concentrated  near  the  Vaisala  frequency. 

The  Garrett  and  Munk  model  spectrum  of  internal  wave  fluctuations  (Garrett  and  Munk, 
1975)  implies  that  the  internal  wave  energy  falls  off  as  o>~7  in  the  frequency  domain.  With  this 
spectrum,  and  for  reasonable  values  of  Vand  /,  90%  of  the  energy  due  to  the  internal  waves  is 
in  a  frequency  band  an  order  of  magnitude  below  the  Vaisala  frequency.  In  this  band,  the  vert¬ 
ical  velocities  all  are  less  than  10%  of  the  horizontal  velocities.  Thus,  we  do  not  do  violence  to 
most  of  the  internal  wave  band  of  motions. 


Now,  it  is  expected  that  dw/dx  and  dw/dy  are  negligible  compared  to  vertical  derivatives 
of  horizontal  velocities  for  mesoscale  and  lower  frequency  motions.  For  internal  waves,  the 
theory  gives 


dw/dx  __  kxw 
du/dz  k2u 


a,2-/2 
~  N1-  a,2’ 


(5) 


For  measured  frequency  spectra,  the  numerator  is  two  orders  of  magnitude  less  than  the 
denominator  in  the  lower  frequency  band  in  which  90%  of  the  energy  resides.  This  implies  that 
the  vertical  shear  is  the  dominant  contributor  to  the  vorticity  in  this  band  also.  This  approxi¬ 
mation  breaks  down  for  the  higher  frequency  waves  for  which  the  wave  frequency  is  order  of  a 
third  the  buoyancy  frequency  or  larger,  and,  in  this  band,  measurements  of  vorticity  are  more 
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useful  than  those  of  horizontal  velocity.  In  principle,  the  ratio  £/N  is  a  measure  of  the  strength 
of  these  wave  motions  in  the  absence  of  fine  structure. 

For  turbulent  motions,  either  in  the  near-surface  active  mixed  layer  or  in  the  microscale 
patches  of  active  turbulence  in  the  interior  of  the  ocean,  all  terms  on  the  right  hand  side  of 
equations  1  are  comparable  in  magnitude.  Since  the  vorticity  is  a  fundamental  turbulence  pro¬ 
perty,  the  measurement  of  it  is  more  useful  in  understanding  the  dynamics  than  is  a  measure¬ 
ment  of  horizontal  velocity.  In  addition,  a  direct  measurement  of  vorticity  near  the  ocean  sur 
face  in  a  search  for  the  turbulence  levels,  momentum  fluxes,  and  the  like  in  the  vicinity  of  the 
surface  has  the  advantage  over  a  velocity  measurement  of  being  practically  independent  of  sur¬ 
face  wave  motions.  The  orbital  motions  of  surface  waves  are  very  nearly  irrotational  so  that 
this  instrument  automatically  suppresses  them. 

Finally,  there  is  a  further  practical  advantage  of  a  vorticity  measurement  over  that  of  velo¬ 
city  because  of  the  shape  of  the  wavenumber  spectrum  of  motions.  What  limited  information 
that  is  available  on  both  the  vertical  and  horizontal  structure  of  motions  in  the  ocean  indicates 
that  the  spectra  are  decidedly  red;  that  is,  spectral  estimates  of  the  longer  waves  are  consider¬ 
ably  higher  than  those  of  the  shorter  waves.  The  spectral  rolloff  is  the  order  of  the  square  of 
the  inverse  wavenumber,  so  that  the  slope  spectra  of  velocity  would  be  more  or  less  white. 
Since  the  vorticity  involves  spatial  derivatives  of  velocities,  it  is  expected  that  spectra  of  the 
vorticity  also  will  be  approximately  white.  Given  a  choice  of  which  variable  to' measure,  the 
vorticity  enjoys  a  considerable  advantage  due  to  its  expected  white  spectrum  because  this  results 
in  higher  reliability  by  the  avoidance  of  side-lobe  effects. 

In  the  assumption  in  which  the  last  terms  in  equation  1  are  neglected,  the  vertical  profile 
of  horizontal  velocity  simply  is  an  integral  of  the  vorticity  profile.  The  velocity  profile  obtained 
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from  the  integral  of  a  measured  vorticity  profile  is  a  relative  one  unless  an  accurate  estimate  of 
water  velocity  at  some  depth  over  the  time  of  the  cast  is  available.  For  a  cast  through  the  per¬ 
manent  thermocline  that  takes  the  order  of  an  hour  to  complete,  present  satellite  or  Loran  C 
navigation  systems  might  be  accurate  enough  to  estimate  ship  drift  for  this  purpose.  For  very 
shallow  casts  through  only  the  mixed  layer  and  seasonal  thermocline,  a  more  accurate  system 
such  as  the  future  NAVSTAR  Global  Positioning  System  or  a  short  range  radio  positioning  sys¬ 
tem  would  be  required. 

Sensor  Considerations 

There  are  many  types  of  instrument  configurations  that  can  sense  water  movement.  Non¬ 
mechanical  ones  receive  high  priority  because  of  their  reliability  and  relative  immunity  to  con¬ 
tamination.  Acoustic  devices,  such  as  the  new  Neil  Brown  current  meter  (N.B.I.S.,  1978),  have 
high  resolution,  accuracy,  and  dynamic  range,  and  they  have  lower  thresholds.  However, 
unless  specially  instrumented,  they  suffer  from  ship  motion. 

The  new  device  considered  here  will  be  called  a  vorticity  meter.  Mechanical  devices  have 
previously  been  designed  to  measure  fluid  vorticity  (cf.  Todd,  1949;  Hopkins  &  Sorensen, 
1956;  Bobrik,  1960;  Lokshminarayana,  1964;  Zalay,  1976;  Wigeland  et  al,  1978;  Govinda  Rtuu 
et  al,  1978),  but  these  all  are  of  the  rotor  or  vane  type.  The  large  number  of  references  on  the 
vane-type  meter  is  due  to  many  attempts  to  construct  responsive,  calibrated  instruments.  This 
technique  still  involves  considerable  art  in  constructing  the  devices,  and  they  undoubtedly 
would  experience  difficulties  if  used  for  any  length  of  time  in  seawater.  Nevertheless,  Gongwer 
and  Finkle  (1962)  have  deployed  a  meter  of  this  type  at  sea  to  assess  the  level  of  "turbulence" 
there.  Although  calibration  and  longevity  are  not  discussed  in  their  article,  the  instrument 
apparently  did  survive  the  sea  test  in  working  order.  Kovasznay  (1954)  has  discussed  an 


7 


J.P.  DUGAN 


arrangement  of  hot  wire  probes  which  measures  fluid  vorticity,  but,  even  in  a  hot  film  instead 
of  hot  wire  arrangement,  these  probes  suffer  calibration  drift  and  high  failure  rates  in  seawater. 
An  acoustic  vorticity  meter  has  been  designed  for  a  towed  device  (Gibson  and  Williams,  1978) 
and  has  been  discussed  in  terms  of  ocean  basin  scale  measurements  (Rossby,  1975),  so  only  its 
use  in  the  profiling  configuration  can  be  considered  to  be  entirely  new.  The  instrument  would 
be  deployed  in  the  same  manner  as  a  CTD  and,  in  its  most  efficient  mode,  both  devices  would 
utilize  the  same  pressure  case  and  electronics.  Two  acoustic  devices  below,  or  on  the  side  of, 
the  case  would  measure  the  two  horizontal  components  of  vorticity.  Each  device  would  be  a 
planar  array  of  projector,  reflectors,  and  receiver,  that  would  measure  the  circulation  around  a 
closed  acoustic  path  in  the  water.  This  measurement  is  accomplished  by  measuring  the  time 
difference  that  a  pulse  takes  to  travel  opposite  directions  around  the  path,  and  relating  this 
directly  to  the  circulation  around  the  path.  This  circulation  T  is  defined  as 

r  -  <£  u  di  (6) 

where  u  is  the  velocity  vector  and  dl  is  the  differential  path  length.  By  Stokes  theorem, 

<J>  u  dl  -  /  £  n  dA  (7) 

where  £  -  V  x  u  is  the  vorticity  vector  and  n dA  is  an  element  of  any  surface  bounded  by  the 

closed  acoustic  path.  Assuming  a  planar  array  and  uniform  distribution  of  vorticity  in  the 
enclosed  area,  this  reduces  to 

T  -  [A,  (8) 

where  £  is  the  vorticity  component  perpendicular  to  the  plane  of  the  array  and  A  is  the  total 

enclosed  area  of  the  array.  Thus  the  vorticity  is  given  by 

£  -  r /A,  (9) 

where  T  is  measured  and  A  is  given  by  the  array  configuration. 


The  travel  time  of  an  acoustic  pulse  (or  constant  phase  point  of  a  CW  signal)  is 
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(10) 


cL  +  (L  u  •  d\ 

JO 

where  L  is  the  path  length,  c  is  the  sound  velocity  and  it  has  been  assumed  that  the  variation  in 
sound  speed  around  the  array  is  small.  The  difference  for  propagation  each  way,  then,  is 


— m) 
where  it  also  has  been  assumed  that  cL  »  |  J*  u  ■  dl|.  Note  that  the  formula  above  does 
indicate  that  the  local  value  of  the  sound  speed  may  be  important.  With  this,  and  the  previous 
formula,  the  vorticity  is  given  by 


£ 


c2A7~ 
2  A  ' 


(12) 


Two  such  devices  yield  two  components  of  horizontal  vorticity,  and,  when  referenced  to  an 
internal  compass,  yield  £x  and  £r  Alternatively,  a  single  array  on  a  spinning  device  with  an 
internal  compass  could  provide  the  required  information.  (See  the  Appendix  for  more  details 
in  the  derivation  of  Eq.  12.) 


There  are  several  critical  design  factors  that  must  be  addressed  before  such  a  system 
would  be  considered  to  be  entirely  feasible.  First,  the  size  of  the  array  defines  the  resolving 
power  of  the  instrument.  The  size  would  be  a  trade-off  between  the  required  resolution  and 
the  accuracy  and  repeatability  of  the  clock  and  pulse  generator,  and  preliminary  estimates  of 
this  size  can  be  made.  The  Richardson  number  defined  by 

R,  -  NVidu/dz)7  (13) 

often  is  between  1  and  100  on  scales  of  meters  to  tens  of  meters  in  the  ocean.  Taking  the 

larger  number  as  representative  of  a  minimum  value  of  vorticity  to  design  for,  and  assuming  a 

typical  buoyancy  frequency  of  10~2  rad/sec  in  the  upper  ocean,  a  typical  noise  floor  of 
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seems  reasonable.  (Note  that  this  could  be  an  order  of  magnitude  smaller  in  the  deeper  ocean.) 
Assuming  a  state-of-the-art  timing  accuracy  of  100  picoseconds  (10~'°  sec)  for  oceanographic 
instrumentation  (Neil  Brown  Acoustic  Current  Meter)  and  a  sound  speed  of  1500  m/s,  equa¬ 
tion  12  allows  a  smallest  acoustic  array  cross  section  of  0.1  m2.  This  corresponds  to  a  triangular 
array  having  legs  of  50  cm.  An  order  of  magnitude  increase  in  timing  accuracy  to  10 
picoseconds  would  allow  this  to  be  reduced  to  about  1 5  cm  on  a  side. 

The  second  problem  area  is  the  level  of  contamination  due  to  instrument  motion.  This 
design  is  entirely  immune  to  translational  motions,  but  it  is  susceptible  to  instrument  rotations. 
Rotation  of  the  acoustic  array  about  an  axis  perpendicular  to  its  own  plane  causes  a  signal  that 
appears  identical  to  fluid  rotation  due  to  vorticity  in  the  form  £  -  2  dO/dt,  Where  d0/dt  is  the 
roll  rate.  This  could  be  subtracted  out  by  proper  instrumentation  such  as  a  roll  rate  gyro  since 
10-3  rad/sec  is  within  that  technology,  but  it  incurs  the  penalty  of  extra  cost  and  complexity. 
This  motion  is  due  to  ship  motion  and  to  the  effects  of  the  shear  on  the  instrument  casing,  and 
it  can  be  modelled  or  tested  experimentally  in  order  to  accurately  assess  its  impact  on  the 
design.  The  crucial  point  is  to  optimize  the  design  of  the  instrument  so  that  the  tilt  effects  are 
minimized. 

A  third  problem  is  the  calculation  of  a  velocity  profile  from  the  vorticity  profile.  For 
many  uses,  this  calculation  is  not  necessary  because  it  is  either  the  vorticity  or  the  vertical  shear 
that  actually  is  required.  For  example,  the  profile  of  Richardson  number,  defined  by  equation 
13  is  an  important  parameter  in  energy  transfer  processes.  However,  for  those  problems  for 
which  the  velocity  is  required,  the  measured  shear  profile  must  be  integrated  in  the  vertical  to 
obtain  the  velocity  profile.  The  resulting  velocity  profile  is  a  relative  one,  and  the  calculation  of 
an  absolute  profile  requires  the  knowledge  of  an  extra  constant— the  absolute  velocity  at  one 
point  in  depth.  This  value  is  not  known  in  general  unless  the  ship  has  an  accurate  navigation 
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system  and  a  2-dimensional  underwater  log.  The  assignment  of  the  resulting  water  velocity  at 
the  depth  of  the  log  ties  down  the  velocity  profile.  There  remains  some  uncertainty  in  the 
accuracy  of  the  velocity  profile  at  depth  since  errors  in  the  measurement  of  shear  could  creep 
up  additively  in  the  integration  procedure. 

Conclusion 

We  have  outlined  a  new  concept  for  measuring  the  Richardson  number  profile  and  the 
velocity  profile  from  survey  vessels  and  have  identified  the  potential  problems  in  its  implemen¬ 
tation.  The  concept  of  measuring  fluid  vorticity  appears  to  be  feasible,  although  solid  estimates 
of  instrument  and  platform  related  noise  sources  are  not  yet  available.  The  proof  of  concept 
rests  on  a  follow  up  engineering  design  study  that  would  emphasize  the  acquisition  of  hard 
numbers  for  the  noise  sources  so  as  to  optimize  the  final  design  options. 
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Appendix 

ASSUMPTIONS  IN  DERIVATION  OF  EQUATION  (12) 

The  equation  governing  the  relationship  between  fluid  vorticity  and  the  sensor  parameters 
is  derived  on  pgs.  8-9.  This  appendix  gives  more  detail  on  the  assumptions  implied  in  that 
derivation. 

The  travel  time  of  an  acoustic  pulse  or  a  constant  phase  point  of  a  CW  signal  over  a  path 
of  length  L  is 


(A.I) 


S'  (c  +  u)  •  d\ 

where  c  is  the  sound  speed,  u  is  the  fluid  velocity  vector,  and  I  is  a  unit  vector  tangent  to  the 
path.  Assuming  that  the  sound  speed  variation  around  the  array  is  small  (but  possibly  the 
sound  speed  is  unknown),  the  equation  simplifies  to 


cL  +  u  ■  d\ 
Jo 

Denoting  the  travel  time  one  way  around  the  path  as  Tu 


(A. 2) 


cL+So*d  1 

the  travel  time  T7  the  other  way  around  the  path  becomes 


(A. 3) 


(A. 4) 


Jo 

The  change  in  the  sign  of  the  integral  in  the  denominator  occurs  because  of  the  change  in 
direction  of  the  vector  I.  The  difference  in  these  times  is 
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and  the  sum  is 


A  T  =  t2-  r, 


2L2r 

c2l2-  r2 


(A. 5) 


T  =  T2  +  T\  =  ■  (A.6) 

where  r  =  f  u  •  </l.  Inspecting  these  two  equations,  it  follows  immediately  that  Eq.  (11) 
results  from  Eq.  (A. 5)  with  the  assumption  that  T  «  cL.  Also,  under  the  same  assumption 
that  T  «  cL ,  it  is  clear  that  the  sum  of  the  travel  times  is  simply  related  to  the  average  sound 
speed  on  the  acoustic  path.  Thus,  assuming  that  both  T  and  the  sound  speed  are  unknown,  and 
that  A  T  and  T  are  measured  quantities,  Eq.  (A. 5)  and  (A. 6)  constitute  two  nonlinear  equations 
for  two  unknowns.  These  equations  have  a  symmetry  which  permits  the  simple  solution 


r  = 


2L2  AT 


cL  =  — 


T2-  AT2 

2 L2  f 


(A. 7) 
(A. 8) 


f 2  -  A  T2 ' 

Thus,  both  c  and  T  can  be  determined  explicitly  from  the  sum  and  difference  of  travel  times. 
The  ratio  AT/T\s  the  order  of  |  u  l/c  which  is  less  than  0.01  for  realistic  values  in  the  ocean, 
so  that 


7"2  -  AT2  -  T2  (1  -  A T2/ T2) 

=  f2 

to  a  very  good  approximation.  In  this  case,  Eqs.  (A. 7)  and  (A. 8)  reduce  to  the  simpler  forms 

r  =  2Z,2A7Vf2  (A. 9) 

c  =  2  LIT.  (A. 10) 

It  is  worth  noting  that  the  assumption  that  the  sound  speed  varies  only  a  small  amount  along 

the  acoustic  path  is  not  necessary,  as  carrying  the  term  cdl  instead  of  simplifying  it  to  cL 

throughout  the  analysis  does  not  change  the  basic  results.  In  that  case,  the  only  change  in  the 

results  is  that  Eq.  (A.  10)  is  replaced  by 
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foL  cdl  -  2LVT.  (A.  11) 

Equation  (12)  follows  from  Eq.  (A.9)  upon  substitution  of  expression  (8)  in  the  latter, 
and  this  exercise  has  underscored  the  accuracy  of  the  assumptions  which  are  implied  in  Eq. 
(12).  This  derivation  also  has  made  it  clear  that  the  instrument  does  not  require  a  separate 
measurement  of  sound  speed  as  a  simple  application  of  Eq.  (12)  would  imply. 


